(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     74442,       2041]
NotebookOptionsPosition[     67996,       1824]
NotebookOutlinePosition[     68396,       1841]
CellTagsIndexPosition[     68353,       1838]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"ClebschGordan", "::", "\"\<tri\>\""}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"ClebschGordan", "::", "\"\<phy\>\""}], "]"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Off", "[", 
   RowBox[{"SixJSymbol", "::", "\"\<tri\>\""}], "]"}], ";"}]}], "Input",
 CellChangeTimes->{{3.456663504453125*^9, 3.45666350465625*^9}, {
  3.5418510742144785`*^9, 3.5418510745114956`*^9}}],

Cell[CellGroupData[{

Cell["Physical Constants", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[
 RowBox[{
  RowBox[{
  "oFile", "=", 
   "\"\<C:\\\\Users\\\\Jason \
Stalnaker\\\\Documents\\\\MathematicaStuff\\\\Rb2Photon\\\\CSEB\\\\\
PowerOpticalDepth\\\\85RbCSEBPOD_30.txt\>\""}], ";", " ", 
  RowBox[{"(*", 
   RowBox[{"output", " ", "file"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.5418532164850097`*^9, 3.5418532218513165`*^9}, {
   3.542081732923125*^9, 3.542081734233525*^9}, {3.545654663150125*^9, 
   3.545654679587625*^9}, {3.545654769775125*^9, 3.54565477097825*^9}, {
   3.54568532929075*^9, 3.54568533035325*^9}, 3.54611592162434*^9, {
   3.54611647515077*^9, 3.546116479782035*^9}, 3.5461165196033125`*^9, {
   3.546119416835025*^9, 3.5461194170780387`*^9}, 3.546123816936696*^9, 
   3.546161262548339*^9, 3.546179174557089*^9, 3.5462541775148044`*^9, 
   3.5462911046825647`*^9, 3.5463360113400784`*^9, {3.5468975716028767`*^9, 
   3.546897584379299*^9}}],

Cell[CellGroupData[{

Cell["Speed of Light", "Subsubtitle"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"c", "=", "29979245800"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"cm", "/", "sec"}], "*)"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Planck's Constant", "Subsubtitle"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"hBar", "=", 
    RowBox[{"1.054", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "27"}]]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"ergs", " ", "sec"}], "*)"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Boltzmann Constant", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"kB", "=", 
    RowBox[{"1.380658", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "16"}]]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"erg", "/", "K"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.4566596001875*^9, 3.456659625953125*^9}, {
  3.45665968075*^9, 3.456659686359375*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Atomic Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["General Parameters", "Subtitle",
 CellChangeTimes->{{3.45666146840625*^9, 3.45666147296875*^9}}],

Cell["Nuclear Spin", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"i", "=", 
   FractionBox["5", "2"]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.4730885615*^9, 3.473088561953125*^9}, 3.53632801824225*^9, 
   3.5386062447474585`*^9, 3.5461791923121047`*^9}],

Cell[CellGroupData[{

Cell["Transit Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.45666389359375*^9, 3.456663897328125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[Gamma]T", "=", 
   FractionBox["1", 
    RowBox[{"60.", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "9"}]]}]]}], ";"}]], "Input",
 CellChangeTimes->{{3.456663726421875*^9, 3.4566637388125*^9}, {
   3.45666380203125*^9, 3.456663813921875*^9}, 3.4566922269375*^9, 
   3.5366835360235*^9, 3.538523997111776*^9, {3.539985422708371*^9, 
   3.5399854227239714`*^9}, {3.54031967759375*^9, 3.5403196779375*^9}, {
   3.541854688631212*^9, 3.541854691954034*^9}, 3.541908707424239*^9, 
   3.5420835944399247`*^9, 3.546161425317891*^9, 3.5468989615159187`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Polarization", "Subsubtitle",
 CellChangeTimes->{{3.456663909203125*^9, 3.456663910671875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"q1", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456663912765625*^9, 3.456663917453125*^9}, 
   3.456663957921875*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"q2", "=", 
   RowBox[{"{", 
    RowBox[{"0", ",", "1", ",", "0"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{3.4566639629375*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Parameters", "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659860875*^9}}],

Cell["Ground State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
   3.456659752203125*^9}, 3.45665980753125*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SGnd", "=", 
   FractionBox["1", "2"]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665975521875*^9, 3.456659757859375*^9}, 3.4566598095625*^9}],

Cell[CellGroupData[{

Cell["Ground State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jGnd", "=", 
   FractionBox["1", "2"]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665975521875*^9, 3.456659757859375*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.456662908015625*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinGnd", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jGnd", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxGnd", "=", 
   RowBox[{"jGnd", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Hyperfine A", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
  3.45666040421875*^9}, {3.456660756828125*^9, 3.45666075759375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"aGnd", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"1011910813.0", ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{"From", " ", "Steck", " ", "review"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.4566605719375*^9}, {3.456660680015625*^9, 
   3.456660705609375*^9}, {3.45666076115625*^9, 3.4566607615625*^9}, {
   3.47308857578125*^9, 3.473088584546875*^9}, {3.47309640846875*^9, 
   3.473096424140625*^9}, {3.536335670945375*^9, 3.53633567542975*^9}, {
   3.53998552417695*^9, 3.539985574127039*^9}, 3.5461792119452276`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ground State Hyperfine B", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
   3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
   3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
   3.45666040421875*^9}, 3.4566605775*^9, {3.45666076390625*^9, 
   3.45666076546875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"bGnd", "=", "0"}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
  3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
  3.4566605581875*^9, 3.456660602203125*^9}, {3.45666063840625*^9, 
  3.45666063971875*^9}, {3.45666076803125*^9, 3.456660773828125*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Parameters", "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659887765625*^9}}],

Cell["Intermediate State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}, {3.45665980753125*^9, 3.45665981840625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SInt", "=", 
   FractionBox["1", "2"]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665975521875*^9, 3.456659757859375*^9}, {3.4566598095625*^9, 
  3.45665982178125*^9}}],

Cell[CellGroupData[{

Cell["Intermediate State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jInt", "=", 
   FractionBox["3", "2"]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.473088588953125*^9, 
   3.47308858946875*^9}, 3.536333582476625*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.456662923203125*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinInt", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jInt", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxInt", "=", 
   RowBox[{"jInt", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
  3.456662927125*^9, 3.456662936640625*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Gamma]Int", "=", 
    FractionBox["1", 
     RowBox[{"26.2", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "9"}]]}]]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Seconds", " ", "From", " ", "Budker"}], ",", "Kimball", ",", " ",
     "DeMille"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.456660078296875*^9, 3.456660092390625*^9}, {3.53668310449225*^9, 
  3.536683106007875*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Center of Gravity Energy", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Nu]Int", "=", 
   RowBox[{"SetPrecision", "[", 
    RowBox[{"384230406373000", ",", "20"}], "]"}]}], " ", 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "from", " ", "Steck"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.473088596515625*^9, 3.473088616875*^9}, {
   3.47308883840625*^9, 3.473088843484375*^9}, {3.473096433640625*^9, 
   3.473096439859375*^9}, {3.53633478836725*^9, 3.53633479311725*^9}, 
   3.536607319211*^9, {3.536607447226625*^9, 3.536607498007875*^9}, {
   3.53666646042975*^9, 3.5366664651485*^9}, {3.537893167558856*^9, 
   3.5378931893520937`*^9}, 3.5378943549667425`*^9, 3.5378974743522215`*^9, {
   3.5378985672711487`*^9, 3.53789856825395*^9}, {3.5378998431661897`*^9, 
   3.537899869374236*^9}, 3.5379017571978636`*^9, {3.5386062999866185`*^9, 
   3.538606301137684*^9}, {3.5399856188991175`*^9, 3.539985624343527*^9}, 
   3.5461795129864464`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Hyperfine A", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
  3.45666040421875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"aInt", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"25002000", ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "From", " ", "Steck"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.4566605719375*^9}, {3.473088627078125*^9, 
   3.4730886284375*^9}, {3.473096443796875*^9, 3.473096448328125*^9}, {
   3.5363357313985*^9, 3.536335734726625*^9}, {3.536607079039125*^9, 
   3.53660709955475*^9}, {3.536666472757875*^9, 3.536666476257875*^9}, {
   3.536683146882875*^9, 3.5366831523985*^9}, {3.5372022874831247`*^9, 
   3.5372022919049997`*^9}, {3.5393622605083046`*^9, 3.5393622891187553`*^9}, 
   3.5461792285501776`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Intermediate State Hyperfine B", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
   3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
   3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
   3.45666013296875*^9, 3.456660152140625*^9}, {3.456660398109375*^9, 
   3.45666040421875*^9}, 3.4566605775*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"bInt", "=", 
    RowBox[{"SetPrecision", "[", 
     RowBox[{"25790000", ",", "20"}], "]"}]}], " ", ";"}], 
  RowBox[{"(*", 
   RowBox[{"Hz", " ", "From", " ", "Steck"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
   3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
   3.456660078296875*^9, 3.456660092390625*^9}, {3.456660306609375*^9, 
   3.45666036596875*^9}, {3.456660407359375*^9, 3.45666043703125*^9}, {
   3.4566605581875*^9, 3.456660602203125*^9}, {3.45666063840625*^9, 
   3.45666063971875*^9}, 3.473088631296875*^9, {3.473096452546875*^9, 
   3.473096456171875*^9}, {3.536335742211*^9, 3.536335745961*^9}, {
   3.536607112711*^9, 3.536607116664125*^9}, 3.53666648311725*^9, {
   3.5366831642735*^9, 3.53668317030475*^9}, {3.5372023036862497`*^9, 
   3.5372023085924997`*^9}, {3.5393622969655685`*^9, 3.5393623013803763`*^9}, 
   3.5461792351735563`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State Parameters", "Subtitle",
 CellChangeTimes->{{3.456659852765625*^9, 3.456659887765625*^9}, {
  3.456660101109375*^9, 3.4566601018125*^9}}],

Cell["Final State S", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.45665974584375*^9, 
  3.456659752203125*^9}, {3.45665980753125*^9, 3.456659829375*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"SFin", "=", 
   FractionBox["1", "2"]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.45665975521875*^9, 3.456659757859375*^9}, {3.4566598095625*^9, 
   3.45665983434375*^9}, {3.536335431601625*^9, 3.536335431726625*^9}, 
   3.53634248861725*^9}],

Cell[CellGroupData[{

Cell["Final State J", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659783546875*^9, 
  3.456659787671875*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"jFin", "=", 
   FractionBox["5", "2"]}], ";"}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
   3.4566597779375*^9, 3.456659780921875*^9}, {3.456661545015625*^9, 
   3.456661545171875*^9}, 3.4569295484212503`*^9, {3.47308863565625*^9, 
   3.473088640640625*^9}, {3.536335439086*^9, 3.536335439226625*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State F Values", "Subsubtitle",
 CellChangeTimes->{{3.45666290321875*^9, 3.45666294725*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fMinFin", "=", 
   RowBox[{"Abs", "[", 
    RowBox[{"jFin", "-", "i"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fMaxFin", "=", 
   RowBox[{"jFin", "+", "i"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456662861609375*^9, 3.456662890515625*^9}, {
  3.456662927125*^9, 3.456662957578125*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Final State Linewidth", "Subsubtitle",
 CellChangeTimes->{{3.456659575640625*^9, 3.456659597484375*^9}, {
  3.45665970596875*^9, 3.456659711984375*^9}, {3.456659762*^9, 
  3.45665976646875*^9}, {3.456659909703125*^9, 3.45665991996875*^9}, {
  3.456659988796875*^9, 3.45665998953125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Gamma]Fin", "=", 
    FractionBox["1", 
     RowBox[{"241.", " ", 
      SuperscriptBox["10", 
       RowBox[{"-", "9"}]]}]]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"Seconds", " ", "From", " ", "Heavens"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.45665971353125*^9, 3.45665971846875*^9}, {
  3.45665976928125*^9, 3.45665977578125*^9}, {3.45665992590625*^9, 
  3.45665992978125*^9}, {3.456659964359375*^9, 3.45665998015625*^9}, {
  3.4566600681875*^9, 3.45666007421875*^9}, {3.4566616303125*^9, 
  3.456661631765625*^9}, {3.4569295617493753`*^9, 3.4569295620150003`*^9}, {
  3.45951917796875*^9, 3.45951917865625*^9}, {3.536683236757875*^9, 
  3.536683244726625*^9}, {3.538582219701305*^9, 3.5385822378273416`*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Energies, Matrix Elements and Transition Probablities", "Subtitle",
 CellChangeTimes->{{3.456681533796875*^9, 3.456681554234375*^9}}],

Cell["Formula for Calculating Energies with HFS", "Subsubtitle",
 CellChangeTimes->{{3.45666156509375*^9, 3.45666157378125*^9}, {
  3.456662570765625*^9, 3.4566625805625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Kappa]", "[", 
    RowBox[{"f_", ",", "j_"}], "]"}], "=", 
   RowBox[{
    RowBox[{"f", 
     RowBox[{"(", 
      RowBox[{"f", "+", "1"}], ")"}]}], "-", 
    RowBox[{"j", 
     RowBox[{"(", 
      RowBox[{"j", "+", "1"}], ")"}]}], "-", 
    RowBox[{"i", 
     RowBox[{"(", 
      RowBox[{"i", "+", "1"}], ")"}]}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.45666169846875*^9, 3.456661713*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Delta]E", "[", 
    RowBox[{"a_", ",", "b_", ",", "f_", ",", "j_"}], "]"}], ":=", 
   RowBox[{"If", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"j", "\[NotEqual]", " ", 
       FractionBox["1", "2"]}], "&&", 
      RowBox[{"i", "\[NotEqual]", 
       FractionBox["1", "2"]}]}], ",", 
     RowBox[{
      RowBox[{
       FractionBox["1", "2"], "a", " ", 
       RowBox[{"\[Kappa]", "[", 
        RowBox[{"f", ",", "j"}], "]"}]}], "+", 
      RowBox[{"b", " ", 
       FractionBox[
        RowBox[{
         RowBox[{
          FractionBox["3", "4"], 
          SuperscriptBox[
           RowBox[{"\[Kappa]", "[", 
            RowBox[{"f", ",", "j"}], "]"}], "2"]}], "+", 
         RowBox[{
          FractionBox["3", "4"], 
          RowBox[{"\[Kappa]", "[", 
           RowBox[{"f", ",", "j"}], "]"}]}], "-", 
         RowBox[{"i", 
          RowBox[{"(", 
           RowBox[{"i", "+", "1"}], ")"}], "j", 
          RowBox[{"(", 
           RowBox[{"j", "+", "1"}], ")"}]}]}], 
        RowBox[{"2", " ", "i", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "i"}], "-", "1"}], ")"}], "j", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "j"}], "-", "1"}], ")"}]}]]}]}], ",", 
     RowBox[{
      FractionBox["1", "2"], "a", " ", 
      RowBox[{"\[Kappa]", "[", 
       RowBox[{"f", ",", "j"}], "]"}]}]}], "]"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456661820734375*^9, 3.45666183171875*^9}, {
   3.456662393875*^9, 3.456662397109375*^9}, {3.456683422046875*^9, 
   3.456683431546875*^9}, 3.460223934290625*^9}],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"\[Nu]GndToInt", "[", 
    RowBox[{"f_", ",", "fPrime_"}], "]"}], "=", 
   RowBox[{"\[Nu]Int", "-", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aGnd", ",", "bGnd", ",", "f", ",", "jGnd"}], "]"}], "+", 
    RowBox[{"\[Delta]E", "[", 
     RowBox[{"aInt", ",", "bInt", ",", "fPrime", ",", "jInt"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Nu]GndToFin", "[", 
   RowBox[{"f_", ",", "fDPrime_"}], "]"}], ":=", 
  RowBox[{"2", " ", 
   RowBox[{"Switch", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"f", ",", "fDPrime"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"3", ",", "5"}], "}"}], ",", "385285142367000", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"3", ",", "4"}], "}"}], ",", "385285147084900", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"3", ",", "3"}], "}"}], ",", "385285151594800", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"3", ",", "2"}], "}"}], ",", "385285155397600", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"3", ",", "1"}], "}"}], ",", "385285158138000", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2", ",", "4"}], "}"}], ",", "385286664961900", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2", ",", "3"}], "}"}], ",", "385286669461000", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2", ",", "2"}], "}"}], ",", "385286673254300", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2", ",", "1"}], "}"}], ",", "385286675998000", ",", 
     "\[IndentingNewLine]", 
     RowBox[{"{", 
      RowBox[{"2", ",", "0"}], "}"}], ",", "385286677451000"}], 
    "]"}]}]}]}], "Input",
 CellChangeTimes->{{3.45686233334375*^9, 3.4568623374375*^9}, 
   3.4572661415462503`*^9, {3.47308885534375*^9, 3.473088860125*^9}, {
   3.4730964853125*^9, 3.473096517875*^9}, {3.473097064765625*^9, 
   3.47309707559375*^9}, {3.536683396664125*^9, 3.53668342911725*^9}, {
   3.5385822740464134`*^9, 3.5385822762455387`*^9}, {3.538582433139513*^9, 
   3.538582435955674*^9}, {3.53858286930646*^9, 3.538583146448312*^9}, {
   3.538606352099599*^9, 3.5386064861562667`*^9}, 3.5414187216357045`*^9}],

Cell[CellGroupData[{

Cell["Reduced Matrix Elements in J Basis", "Subsubtitle"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"dRedJGndToInt", "=", "1"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{"(", "4.51", ")"}], "1.28", " ", 
     SuperscriptBox["10", "6"]}], " ", 
    RowBox[{"(*", 
     RowBox[{"Hz", "/", 
      RowBox[{"(", 
       RowBox[{"V", "/", "cm"}], ")"}]}], "*)"}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"dRedJIntToFin", "=", "1"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     FractionBox["1", 
      RowBox[{"e", " ", "a0"}]], 
     RowBox[{"\[Sqrt]", 
      RowBox[{"(", 
       RowBox[{
        FractionBox[
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", 
             FractionBox["1", "2"]}], "+", "1"}], " ", ")"}], "3", " ", 
          "hBar", " ", 
          SuperscriptBox["c", "3"]}], 
         RowBox[{"4", " ", 
          SuperscriptBox[
           RowBox[{"(", 
            RowBox[{"2", " ", "\[Pi]", " ", "\[Nu]6P12To8S"}], ")"}], "3"]}]],
         "2.04", " ", 
        SuperscriptBox["10", "6"]}], ")"}]}], "1.28", " ", 
     SuperscriptBox["10", "6"]}], " ", 
    RowBox[{"(*", 
     RowBox[{"Hz", "/", 
      RowBox[{"(", 
       RowBox[{"V", "/", "cm"}], ")"}]}], "*)"}], ";"}], "*)"}]}]}], "Input",
 CellChangeTimes->{{3.4566056364375*^9, 3.4566056744375*^9}, {
  3.456662686453125*^9, 3.456662724453125*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Reduced Matrix Elements in F Basis", "Subsubtitle",
 CellChangeTimes->{3.4566632656875*^9}],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      RowBox[{
       SubscriptBox["dRedGndToInt", 
        RowBox[{"f", ",", "fPrime"}]], "=", 
       RowBox[{"\[Sqrt]", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", " ", "fPrime"}], "+", "1"}], ")"}], 
          RowBox[{"(", 
           RowBox[{
            RowBox[{"2", "f"}], "+", "1"}], ")"}], 
          SuperscriptBox[
           RowBox[{"SixJSymbol", "[", 
            RowBox[{
             RowBox[{"{", 
              RowBox[{"jInt", ",", "fPrime", ",", "i"}], "}"}], ",", 
             RowBox[{"{", 
              RowBox[{"f", ",", "jGnd", ",", "1"}], "}"}]}], "]"}], "2"], 
          SuperscriptBox["dRedJGndToInt", "2"]}], ")"}]}]}], ";"}], ",", 
     RowBox[{"{", 
      RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], ",",
    "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"f", ",", "fMinGnd", ",", "fMaxGnd"}], "}"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.456662795671875*^9, 3.456662835546875*^9}, {
  3.456663002453125*^9, 3.456663020171875*^9}, {3.5372063527331247`*^9, 
  3.5372063614518747`*^9}}],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Do", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{
      SubscriptBox["dRedIntToFin", 
       RowBox[{"fPrime", ",", "fDPrime"}]], "=", 
      RowBox[{"\[Sqrt]", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", " ", "fDPrime"}], "+", "1"}], ")"}], 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"2", "fPrime"}], "+", "1"}], ")"}], 
         SuperscriptBox[
          RowBox[{"SixJSymbol", "[", 
           RowBox[{
            RowBox[{"{", 
             RowBox[{"jFin", ",", "fDPrime", ",", "i"}], "}"}], ",", 
            RowBox[{"{", 
             RowBox[{"fPrime", ",", "jInt", ",", "1"}], "}"}]}], "]"}], "2"], 
         SuperscriptBox["dRedJIntToFin", "2"]}], ")"}]}]}], ",", 
     RowBox[{"{", 
      RowBox[{"fDPrime", ",", "fMinFin", ",", "fMaxFin"}], "}"}]}], "]"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.456663051265625*^9, 3.456663093375*^9}, {
   3.456663125046875*^9, 3.456663154296875*^9}, 3.456663192609375*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Matrix Elements in F Basis", "Subsubtitle",
 CellChangeTimes->{{3.456663276328125*^9, 3.456663278140625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Do", "[", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Do", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           SubscriptBox["dGndToInt", 
            RowBox[{"f", ",", "fPrime", ",", "m", ",", "mPrime"}]], "=", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], 
             RowBox[{"fPrime", "-", "mPrime"}]], 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"n", "=", "1"}], "3"], 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                "q1", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
                RowBox[{"ThreeJSymbol", "[", 
                 RowBox[{
                  RowBox[{"{", 
                   RowBox[{"fPrime", ",", 
                    RowBox[{"-", "mPrime"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"1", ",", 
                    RowBox[{"n", "-", "2"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"f", ",", "m"}], "}"}]}], "]"}]}], ")"}], 
              SubscriptBox["dRedGndToInt", 
               RowBox[{"f", ",", "fPrime"}]]}]}]}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{"m", ",", 
            RowBox[{"-", "f"}], ",", "f"}], "}"}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"mPrime", ",", 
          RowBox[{"-", "fPrime"}], ",", "fPrime"}], "}"}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"f", ",", "fMinGnd", ",", "fMaxGnd"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.45666329196875*^9, 3.4566633891875*^9}, {3.456663993*^9,
   3.45666408128125*^9}, {3.45666411909375*^9, 3.456664121515625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Do", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Do", "[", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"Do", "[", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"Do", "[", "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{
           SubscriptBox["dIntToFin", 
            RowBox[{
            "fPrime", ",", "fDPrime", ",", "mPrime", ",", "mDPrime"}]], "=", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{"-", "1"}], ")"}], 
             RowBox[{"fDPrime", "-", "mDPrime"}]], 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"n", "=", "1"}], "3"], 
             RowBox[{
              RowBox[{"(", 
               RowBox[{
                RowBox[{
                "q2", "\[LeftDoubleBracket]", "n", "\[RightDoubleBracket]"}], 
                RowBox[{"ThreeJSymbol", "[", 
                 RowBox[{
                  RowBox[{"{", 
                   RowBox[{"fDPrime", ",", 
                    RowBox[{"-", "mDPrime"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"1", ",", 
                    RowBox[{"n", "-", "2"}]}], "}"}], ",", 
                  RowBox[{"{", 
                   RowBox[{"fPrime", ",", "mPrime"}], "}"}]}], "]"}]}], ")"}], 
              SubscriptBox["dRedIntToFin", 
               RowBox[{"fPrime", ",", "fDPrime"}]]}]}]}]}], ",", 
          "\[IndentingNewLine]", 
          RowBox[{"{", 
           RowBox[{"mPrime", ",", 
            RowBox[{"-", "fPrime"}], ",", "fPrime"}], "}"}]}], "]"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"{", 
         RowBox[{"mDPrime", ",", 
          RowBox[{"-", "fDPrime"}], ",", "fDPrime"}], "}"}]}], "]"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{"{", 
       RowBox[{"fDPrime", ",", "fMinFin", ",", "fMaxFin"}], "}"}]}], "]"}], 
    ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.45666329196875*^9, 3.456663489203125*^9}, {
  3.456664105703125*^9, 3.45666415221875*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Doppler Distribution Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"mRb", "=", 
   RowBox[{"1.41", " ", 
    SuperscriptBox["10", 
     RowBox[{"-", "22"}]]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"temp", "=", "323"}], ";"}]}], "Input",
 CellChangeTimes->{{3.456606503203125*^9, 3.4566065051875*^9}, {
   3.456664763890625*^9, 3.456664790515625*^9}, {3.473097352640625*^9, 
   3.473097368078125*^9}, {3.536605457211*^9, 3.53660545749225*^9}, 
   3.5366060920235*^9, {3.536606131101625*^9, 3.536606151414125*^9}, 
   3.53666517880475*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[Alpha]", "=", 
  RowBox[{"\[Sqrt]", 
   RowBox[{"(", 
    FractionBox[
     RowBox[{"2", " ", "kB", " ", "temp"}], "mRb"], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.45666479665625*^9, 3.456664796984375*^9}, {
  3.536606154070375*^9, 3.536606154336*^9}}],

Cell[BoxData["25150.681617316117`"], "Output",
 CellChangeTimes->{3.5420839275623245`*^9, 3.5424503045508747`*^9, 
  3.545654731431375*^9, 3.545654785368875*^9, 3.545685374900125*^9, 
  3.54611605257097*^9, 3.5461207561036267`*^9, 3.5461238608922105`*^9, 
  3.5461615405954843`*^9, 3.5461795733779*^9, 3.546254222444374*^9, 
  3.54629113275117*^9, 3.5463360432649045`*^9, 3.5468976141909513`*^9, 
  3.5468989781465483`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"dopplerDist", "[", "v_", "]"}], ":=", 
  RowBox[{"Exp", "[", 
   RowBox[{"-", 
    FractionBox[
     SuperscriptBox["v", "2"], 
     SuperscriptBox["\[Alpha]", "2"]]}], "]"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Parameters", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell["Offset Frequency", "Subsubtitle"],

Cell[BoxData[
 RowBox[{
  RowBox[{"f0", "=", 
   RowBox[{
    RowBox[{"-", "40."}], " ", 
    SuperscriptBox["10", "6"]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.536606229336*^9, 3.536606229476625*^9}, 
   3.536609033945375*^9, {3.53667308924225*^9, 3.536673089507875*^9}, {
   3.536679694632875*^9, 3.5366796956485*^9}, {3.536683666961*^9, 
   3.5366836672735*^9}, 3.5370104677487497`*^9, 3.5370325434049997`*^9, {
   3.539116492685813*^9, 3.5391164930290136`*^9}, 3.539127919839781*^9, {
   3.540137940370176*^9, 3.540137940652192*^9}, {3.540230418171875*^9, 
   3.54023042025*^9}, {3.545654603587625*^9, 3.545654603900125*^9}, {
   3.546897618371759*^9, 3.54689761888656*^9}}],

Cell[CellGroupData[{

Cell["Rep Rate Range", "Subsubtitle",
 CellChangeTimes->{{3.45666493915625*^9, 3.456664941796875*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"fRepMin", "=", 
   RowBox[{"250.", " ", 
    SuperscriptBox["10", "6"]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fRepMax", "=", 
   RowBox[{
    RowBox[{"250.", " ", 
     SuperscriptBox["10", "6"]}], "+", "90"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.456664947171875*^9, 3.4566649499375*^9}, {
   3.456665042671875*^9, 3.45666507275*^9}, {3.456665103265625*^9, 
   3.45666510571875*^9}, {3.456665200171875*^9, 3.45666520296875*^9}, {
   3.456691746078125*^9, 3.45669175240625*^9}, {3.536342589414125*^9, 
   3.536342618445375*^9}, {3.5366062423985*^9, 3.536606269101625*^9}, {
   3.536673100086*^9, 3.536673151445375*^9}, {3.536679707945375*^9, 
   3.536679725882875*^9}, {3.536683688101625*^9, 3.53668372074225*^9}, {
   3.5369246547956247`*^9, 3.5369246550456247`*^9}, {3.5372003428581247`*^9, 
   3.5372003561081247`*^9}, {3.5377186875620003`*^9, 
   3.5377186920307503`*^9}, {3.5377228732026253`*^9, 
   3.5377228764838753`*^9}, {3.5377945896870003`*^9, 
   3.5377945969526253`*^9}, {3.5378884679682503`*^9, 
   3.5378884728745003`*^9}, {3.5378886446870003`*^9, 
   3.5378886463745003`*^9}, {3.5378933335587473`*^9, 
   3.5378933356803513`*^9}, {3.537893370718013*^9, 3.5378933755228214`*^9}, {
   3.5379017917519245`*^9, 3.537901792812726*^9}, {3.5385833273276577`*^9, 
   3.538583332126932*^9}, {3.539116509315442*^9, 3.539116540281497*^9}, 
   3.5401379532159104`*^9, {3.5401550054482765`*^9, 3.5401550066962786`*^9}, {
   3.540155045134746*^9, 3.540155077052402*^9}, {3.540230423015625*^9, 
   3.54023043365625*^9}, {3.54032197040625*^9, 3.540321972640625*^9}, {
   3.541854003677622*^9, 3.541854008856855*^9}, {3.541907018817148*^9, 
   3.5419070190199547`*^9}, {3.5420814821375246`*^9, 3.542081482808325*^9}, {
   3.54565442847825*^9, 3.54565444879075*^9}, {3.54565449697825*^9, 
   3.5456545043845*^9}, {3.54568535972825*^9, 3.54568536885325*^9}, {
   3.5461162315408363`*^9, 3.546116235548066*^9}, {3.5461163826674805`*^9, 
   3.546116397237314*^9}, {3.546161329652419*^9, 3.546161333378632*^9}, {
   3.546897635484989*^9, 3.5468976432694025`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"fRepAvg", "=", 
   FractionBox[
    RowBox[{"fRepMin", "+", "fRepMax"}], "2"]}], ";"}]], "Input",
 CellChangeTimes->{{3.456665192578125*^9, 3.456665212796875*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Modes", "Subsubtitle",
 CellChangeTimes->{{3.4566813895*^9, 3.45668139175*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n1NomA", "=", 
   RowBox[{"Round", "[", 
    FractionBox[
     RowBox[{"\[Nu]GndToInt", "[", 
      RowBox[{"fMinGnd", ",", "fMinInt"}], "]"}], "fRepAvg"], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n1NomB", "=", 
   RowBox[{"Round", "[", 
    FractionBox[
     RowBox[{"\[Nu]GndToInt", "[", 
      RowBox[{"fMaxGnd", ",", "fMaxInt"}], "]"}], "fRepAvg"], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.456665147765625*^9, 3.456665238296875*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"n2NomA", "=", 
   RowBox[{"Round", "[", 
    RowBox[{
     FractionBox["1", "fRepAvg"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Nu]GndToFin", "[", 
        RowBox[{"fMinGnd", ",", "fMinFin"}], "]"}], "-", 
       RowBox[{"\[Nu]GndToInt", "[", 
        RowBox[{"fMinGnd", ",", "fMinInt"}], "]"}]}], ")"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"n2NomB", "=", 
   RowBox[{"Round", "[", 
    RowBox[{
     FractionBox["1", "fRepAvg"], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"\[Nu]GndToFin", "[", 
        RowBox[{"fMaxGnd", ",", "fMaxFin"}], "]"}], "-", 
       RowBox[{"\[Nu]GndToInt", "[", 
        RowBox[{"fMaxGnd", ",", "fMaxInt"}], "]"}]}], ")"}]}], "]"}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.456665147765625*^9, 3.4566653128125*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n1List", "=", 
  RowBox[{"Union", "[", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n1NomA", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n1NomB", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}]}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.456665383359375*^9, 3.45666546196875*^9}, {
  3.456665618125*^9, 3.456665620609375*^9}, {3.5372066701081247`*^9, 
  3.5372066730768747`*^9}, {3.5377187815151253`*^9, 3.5377187956088753`*^9}, {
  3.5377946040932503`*^9, 3.5377946077338753`*^9}, {3.5377949385776253`*^9, 
  3.5377949450151253`*^9}, {3.537893207650926*^9, 3.537893217432143*^9}, {
  3.5385234618671618`*^9, 3.53852346744048*^9}, {3.538583211301021*^9, 
  3.5385832213905983`*^9}, {3.540321985453125*^9, 3.540321988328125*^9}, {
  3.546116404195712*^9, 3.546116415913382*^9}, {3.5468976521458187`*^9, 
  3.546897665281041*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1536913", ",", "1536914", ",", "1536915", ",", "1536916", ",", "1536917", 
   ",", "1536918", ",", "1536919", ",", "1536920", ",", "1536921", ",", 
   "1536924", ",", "1536925", ",", "1536926", ",", "1536927", ",", "1536928", 
   ",", "1536929", ",", "1536930", ",", "1536931", ",", "1536932"}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.542083927609125*^9, 3.5424503045821247`*^9, 3.54565473147825*^9, 
   3.54565478541575*^9, 3.545685374962625*^9, 3.54611605261777*^9, 
   3.5461162403453403`*^9, {3.5461164015275593`*^9, 3.546116416349407*^9}, 
   3.546120756145629*^9, 3.5461238609152117`*^9, 3.546161540627486*^9, 
   3.5461795775191374`*^9, 3.546254222477376*^9, 3.5462911327921724`*^9, 
   3.546336043293906*^9, 3.5468976664198437`*^9, 3.5468989781777477`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"n2List", "=", 
  RowBox[{"Union", "[", 
   RowBox[{"Join", "[", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n2NomA", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}], ",", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"n2NomB", "+", "\[Delta]n"}], ",", 
       RowBox[{"{", 
        RowBox[{"\[Delta]n", ",", 
         RowBox[{"-", "4"}], ",", "4"}], "}"}]}], "]"}]}], "]"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.456665383359375*^9, 3.45666546196875*^9}, {
  3.456665585125*^9, 3.45666558903125*^9}, {3.456665624328125*^9, 
  3.456665626734375*^9}, {3.5372066753268747`*^9, 3.5372066804518747`*^9}, {
  3.5377187866401253`*^9, 3.5377187978745003`*^9}, {3.5377946088901253`*^9, 
  3.5377946116245003`*^9}, {3.5377949469370003`*^9, 3.5377949552338753`*^9}, {
  3.5378932198813477`*^9, 3.537893225029357*^9}, {3.538523470201638*^9, 
  3.5385234737538414`*^9}, {3.538583224257762*^9, 3.5385832272179317`*^9}, {
  3.540321990921875*^9, 3.54032199415625*^9}, {3.5461164084199533`*^9, 
  3.546116422729772*^9}, {3.5468976686974473`*^9, 3.546897672020253*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "1545360", ",", "1545361", ",", "1545362", ",", "1545363", ",", "1545364", 
   ",", "1545365", ",", "1545366", ",", "1545367", ",", "1545368", ",", 
   "1545369"}], "}"}]], "Output",
 CellChangeTimes->{3.542083927624725*^9, 3.5424503045977497`*^9, 
  3.545654731493875*^9, 3.54565478541575*^9, 3.54568537497825*^9, 
  3.5461160526333704`*^9, 3.5461162423224535`*^9, 3.546116423029789*^9, 
  3.54612075616063*^9, 3.5461238609272127`*^9, 3.5461615406404867`*^9, 
  3.546179577535138*^9, 3.546254222491377*^9, 3.546291132805173*^9, 
  3.5463360433089075`*^9, 3.5468976723946543`*^9, 3.546898978193348*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comb Equation", "Subsubtitle",
 CellChangeTimes->{{3.456681410328125*^9, 3.456681413015625*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Omega]L", "[", 
    RowBox[{"n_", ",", "fR_"}], "]"}], ":=", 
   RowBox[{"2", " ", "\[Pi]", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"n", " ", "fR"}], "+", "f0"}], ")"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.456678515171875*^9, 3.45667856759375*^9}, {
  3.4566795630625*^9, 3.45667956609375*^9}}],

Cell[CellGroupData[{

Cell["Power Optical Depth", "Section",
 CellChangeTimes->{{3.5456539738845*^9, 3.54565397597825*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"NumberDensity", " ", "=", " ", 
    SuperscriptBox["10", "11"]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{"atoms", " ", "per", " ", "cubic", " ", "cm"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.54547662804075*^9, 3.545476669150125*^9}, {
  3.545646682431375*^9, 3.545646692212625*^9}, {3.545647538900125*^9, 
  3.5456475391345*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Lambda]", "=", 
    RowBox[{"780.", " ", 
     SuperscriptBox["10", 
      RowBox[{"-", "7"}]]}]}], ";"}], " ", 
  RowBox[{"(*", 
   RowBox[{
   "Wavelength", " ", "for", " ", "D2", " ", "Transition", " ", "in", " ", 
    "cm"}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.546703806379816*^9, 3.5467038204666414`*^9}, 
   3.546703918294413*^9, 3.5467039561400795`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"\[Delta]v", "=", 
    RowBox[{
     FractionBox[
      RowBox[{"\[Gamma]Int", "+", "\[Gamma]T"}], 
      RowBox[{"2", " ", "\[Pi]"}]], "\[Lambda]"}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"Natural", " ", "Linewidth", " ", "in", " ", 
    RowBox[{"cm", "/", "s"}]}], "*)"}]}]], "Input",
 CellChangeTimes->{{3.54564749829075*^9, 3.545647520572*^9}, {
  3.546602651970539*^9, 3.546602655106145*^9}, {3.54660324445958*^9, 
  3.546603246066383*^9}, {3.546603279044841*^9, 3.546603285784053*^9}, {
  3.546603359494182*^9, 3.546603382083022*^9}, {3.546603430037506*^9, 
  3.5466034342807136`*^9}, {3.546703794508196*^9, 3.5467038041022124`*^9}, {
  3.5467041253223767`*^9, 3.546704126336378*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Table", "[", 
  RowBox[{
   RowBox[{
    SubscriptBox["\[Sigma]", 
     RowBox[{"f_", ",", "fPrime_"}]], "=", 
    RowBox[{
     FractionBox[
      SuperscriptBox["\[Lambda]", "2"], 
      RowBox[{"2", " ", "\[Pi]"}]], 
     RowBox[{"(", 
      FractionBox[
       RowBox[{
        RowBox[{"2", " ", "fPrime"}], "+", "1"}], 
       RowBox[{
        RowBox[{"2", "f"}], "+", "1"}]], ")"}], 
     FractionBox[
      SuperscriptBox[
       SubscriptBox["dRedGndToInt", 
        RowBox[{"f", ",", "fPrime"}]], "2"], 
      RowBox[{
       SuperscriptBox[
        SubscriptBox["dRedGndToInt", 
         RowBox[{"fMaxGnd", ",", "fPrime"}]], "2"], "+", 
       SuperscriptBox[
        SubscriptBox["dRedGndToInt", 
         RowBox[{"fMinGnd", ",", "fPrime"}]], "2"]}]]}]}], ",", 
   RowBox[{"{", 
    RowBox[{"f", ",", "fMinGnd", ",", "fMaxGnd"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"fPrime", ",", "fMinInt", ",", "fMaxInt"}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.546702951513915*^9, 3.546703055831298*^9}, 
   3.5467030890593567`*^9, 3.5467031695710983`*^9, {3.546703256510051*^9, 
   3.5467033302201805`*^9}, {3.5467033647742405`*^9, 3.546703383104273*^9}, {
   3.54670348973046*^9, 3.5467034971872735`*^9}, {3.5467035684325986`*^9, 
   3.5467036389603224`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
    "5.809792042626548`*^-10", ",", "7.531211907108488`*^-10", ",", 
     "6.02496952568679`*^-10", ",", "0"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
    "0", ",", "1.536982021858875`*^-10", ",", "5.379437076506064`*^-10", ",", 
     "1.244955437705689`*^-9"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{
  3.546703882819951*^9, {3.546703938839649*^9, 3.546703958230483*^9}, 
   3.546704127584381*^9, 3.546704811580497*^9, 3.5467090987167068`*^9, 
   3.5467133024772096`*^9, 3.5467133882461605`*^9, 3.546713700869688*^9, 
   3.5467166881195493`*^9, 3.5467183312885327`*^9, 3.5467199297799616`*^9, 
   3.5467201880077314`*^9, 3.5467232506849065`*^9, 3.546725065411703*^9, 
   3.5467270846421957`*^9, 3.546729377349332*^9, 3.5467295439978633`*^9, 
   3.546897749111989*^9, 3.5468989782245483`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"n", "[", "v_", "]"}], "=", 
   RowBox[{
    FractionBox["NumberDensity", 
     RowBox[{"\[Alpha]", " ", 
      RowBox[{"\[Sqrt]", "\[Pi]"}]}]], " ", 
    RowBox[{"dopplerDist", "[", "v", "]"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.545476674650125*^9, 3.54547672816575*^9}, {
  3.545478129322*^9, 3.54547813004075*^9}, {3.5456466452595*^9, 
  3.5456466502595*^9}, {3.545648680806375*^9, 3.54564868172825*^9}, {
  3.5467030936301646`*^9, 3.5467031180286074`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"NRes", "[", "v_", "]"}], ":=", 
  RowBox[{"Integrate", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"n", "[", "vPrime", "]"}], 
     FractionBox[
      FractionBox[
       SuperscriptBox["\[Delta]v", "2"], "4"], 
      RowBox[{
       FractionBox[
        SuperscriptBox["\[Delta]v", "2"], "4"], "+", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"v", "-", "vPrime"}], ")"}], "2"]}]]}], ",", 
    RowBox[{"{", 
     RowBox[{"vPrime", ",", 
      RowBox[{"-", "\[Infinity]"}], ",", "\[Infinity]"}], "}"}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.545648686118875*^9, 3.54564875841575*^9}, 
   3.545648790743875*^9, {3.5456489971345*^9, 3.545649023712625*^9}, {
   3.545649061806375*^9, 3.5456490620095*^9}, 3.545649278025125*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Timing", "[", 
  RowBox[{
   RowBox[{"voigtList", "=", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"v", ",", 
        RowBox[{"NRes", "[", "v", "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"v", ",", 
        RowBox[{"-", "100000"}], ",", "100000", ",", "1000"}], "}"}]}], 
     "]"}]}], ";"}], "]"}]], "Input"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"113.179`", ",", "Null"}], "}"}]], "Output",
 CellChangeTimes->{3.545654761900125*^9, 3.545654924368875*^9, 
  3.54568549291575*^9, 3.5461161780767784`*^9, 3.546120871187209*^9, 
  3.5461239724745927`*^9, 3.5461616592062683`*^9, 3.5461796949568543`*^9, 
  3.5462543389582233`*^9, 3.5462912648297243`*^9, 3.546336167373003*^9, 
  3.5466032054283113`*^9, 3.546603535540491*^9, 3.54660481531714*^9, 
  3.5466096630100546`*^9, 3.546613044393994*^9, 3.54661664626872*^9, 
  3.546622235747343*^9, 3.546703904269988*^9, 3.546703953472475*^9, 
  3.5467040932175207`*^9, 3.546704219858543*^9, 3.546704931466354*^9, 
  3.546709216016416*^9, 3.5467133467656875`*^9, 3.5467138442578897`*^9, 
  3.5467168008159947`*^9, 3.546718442976921*^9, 3.5467200432694526`*^9, 
  3.5467203070045376`*^9, 3.546723379464272*^9, 3.546725186031602*^9, 
  3.5467272187998695`*^9, 3.54672949551309*^9, 3.546729668433981*^9, 
  3.546897868171398*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"NFn", "=", 
   RowBox[{"Interpolation", "[", "voigtList", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"ListPlot", "[", "voigtList", "]"}]}], "Input",
 CellChangeTimes->{
  3.546603106212137*^9, {3.546603210981921*^9, 3.546603212151923*^9}}],

Cell[BoxData[
 GraphicsBox[
  {Hue[0.67, 0.6, 0.6], PointBox[CompressedData["
1:eJxVVwlQ1HUUBkfSJEm8MBWJEo9QC8xjFOVxKl4IaeAFcUnSEKsGom6GONoq
uGoiKuVRjre4OZIi5xMUJE5hYRdQYWEFORb2XrCMrN/PprczOzvffLO7v997
33vf97cPjfGPGGRmZlby+v3P57+vqD4s04xMxakL3P6BCY59OMIsTT2kgGEz
lQl91kkrnEMXMv6mCf1yxjScH+rCeIEJS04df5nzG8MJTibM3jd878CWRYzX
GTHaFPNY7bCY8RlGHPTg1xszOhk2izOiykmyePNdV8bPM+IlkZ947hJgfL8B
U8OUGe9lMpyQZcBhk9fe6pzlxnihAaXZHSv9bjCcsMiAUzD58mgnd8YP6NEi
fr/1B3kMJ6AemyrCYk+s8WB8oh7n1WN8oYHhBA89Wh36KM/rgifjLfQoahNe
eDvYi/FFOtyVbt4T4ujNeJEOhWJ56wTLJYz30WHeFpfQU+ZLGW+pwzqF8WjH
SB/Gl2lxckS40uCxjPFiLb4au/GTn9OWM95Xi1Jv8dk/bVcy3lqLLsMLtuyp
W8X4ag3eyDc42EeuZnyKBs27Z5ReFfkxfq0G69ZZ9JnJ/Blvo0HvR6kN59ev
YbxcjVXFy2KTbT9nfJoaw5p9a5IdAxm/QY1eBWfWx/6wnvG2apRYwaGY0E0M
5/biwq+dZkwN+IJh6EXh9Tnn9V6hDBf24KRX97a5Hgpn2LsHAz60Fh5xjmS4
RIUyr5zIw75RDK9QYavfJtUifTTDld04PPemxG7CVob9u7G/ujp88MntDNd2
oWzsvZqU3XEMB3bhzGNDEzzlOxlu7MTiSSPAL+9bhoM6scV1zsuY9/cyrOjA
VOlZl6SsfQyHd+DxpgmKcfcOMNz+Ao15wbtLpx1iOOoF6r9PcoixEfN5aMct
wZIjBp9jXP/t2LssPM5OdpzrvQ1HXwmdlJudyvXdhrMP6gUyszSu5+c43DNz
TI/kDNfvc/zKcW7eNx3nuV6VGFYpFw5ecIHrU4lTZV4b2youcj0qse2zgSvR
6Ve4/lqx0SN9h6bhOtdbK5Ylu/ZnbZRwfbXg2ovBFX0/3uJ6akFfq3eSCzJv
c/0oMO+PRBw25A7XiwJPJk2Ud4syuT6acdOoybbukM310IxJJ5rSj8/K4/1v
wsP5+dYHriHv9zN0zC1J+XR9Ae/vU3zilnIpbsoD3s8nuLzv1LTp7xbx/jXi
Ns/Emc9GPeL9asDAqNQpstm/8/7UY5HATToQXcb7Icdy6ZcOnVjB6y/D1KuK
kxEnqni96zA2MiQjR/qY17cWNyQF2M+1reH1lKJ1/MeXgqKlvH41eNS5LPZE
YS2vVzWGdK0KsZok4/V5jD1J8UrpLjmvRxUW1VtYVsnq+f0r8U5D7PRi50Z+
33L8SScOijz4hN+vFAMipeU+9U/5fUpw/l82e17aNfHzF+P8rT6K72Y28/M+
xKQZK4dqCjkWFaLPHrsR4/0V/Dz3Mf266HR5HceQj1XXzNO7Vrfw/8tGaYV5
SUQ+x4K7mDZ3abG7fSv/vdsoGTdq/y+xHIMET++0Wng/m2PBZdygmrj9pfYN
fw7fuqBOKh2v5PgI7nW/mVPrzPHr1+ZtOw43LfiPB8LDOSDfF1wG+vsSIP8v
ug3kfIK7QM4fmA3kfpAP5P4p94HUR1QIpH7Ch0DqKygGUv/wEiD9CSwF0r8V
5UD6C5VA+p9WBUQfKY+B6EdcDURfohog+kuUAtGnsBaIfuPqgOhbIAOi/yg5
kPkIrwcyP0ENQOYrsBHI/Pk/ATKfK54CmV/vZ0DmG5qAzL9tM5D9kNYMZH/Y
KIDslxQFkP1j3QJkP4lbgOwvy1Yg+03UCmT/WSiB7MdEJZD9OaAEsl+Fz4Hs
3/7nQPZzXBuQ/a1rA7LfBe1A9r+qHYg/RL0A4h/tL4D4S3gHEP9RdADxp6BO
IP7V2AnE3wK7gPhfbRcQf/TvBuKfld1A/HWFCoj/lqiA+LN3DxD/LuwB4u/Q
C8T/c3uB5ANbNfw/P7zOF0DyRZoaSP6Qq4HkExsNkPyyVgMk36RogOSfag2Q
fGStBZKffLVA8pVYCyR/lWmB5DNLHZD85qMDku9EOiD5r0gHJB9a6IHkRw89
kHyZqAeSP1EPJJ8O6IHk10UGIPlWaACSf7MMQPJxvwFIfp5nBJKv44xA8neG
EUg+1xmB5HcnE5B8LzAByf83TUCeD1QmIM8Pjn1Ani+i+uDN88ffRzmHDQ==

    "]]},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->True,
  PlotRange->Automatic,
  PlotRangeClipping->True]], "Output",
 CellChangeTimes->{3.5466032055063114`*^9, 3.5466035355716915`*^9, 
  3.5466036795131445`*^9, 3.546609663025655*^9, 3.5466130444095936`*^9, 
  3.54661664628432*^9, 3.5466222357941437`*^9, 3.5467040933111205`*^9, 
  3.546704222822548*^9, 3.5467049315633593`*^9, 3.546709216174425*^9, 
  3.546713844431899*^9, 3.546716800841996*^9, 3.5467184430039225`*^9, 
  3.546720043309455*^9, 3.54672030704954*^9, 3.5467233794922733`*^9, 
  3.5467251860666037`*^9, 3.5467272188838744`*^9, 3.5467294955710936`*^9, 
  3.5467296685389867`*^9, 3.5468978682337985`*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"R", "[", 
    RowBox[{"L_", ",", "v_", ",", "f_", ",", "fPrime_"}], "]"}], "=", 
   RowBox[{"Exp", "[", 
    RowBox[{
     RowBox[{"-", 
      SubscriptBox["\[Sigma]", 
       RowBox[{"f", ",", "fPrime"}]]}], " ", "0.28", "  ", 
     FractionBox[
      RowBox[{
       RowBox[{"2", "f"}], "+", "1"}], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", "fMinGnd"}], "+", "1"}], ")"}], "+", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{"2", "fMaxGnd"}], "+", "1"}], ")"}]}]], 
     RowBox[{"NFn", "[", "v", "]"}], " ", "L"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"length", "=", "3.0"}], ";", 
  RowBox[{"(*", 
   RowBox[{"Effective", " ", "Optical", " ", "Depth", " ", "in", " ", "cm"}], 
   "*)"}]}]}], "Input",
 CellChangeTimes->{{3.545653895681375*^9, 3.54565390729075*^9}, {
   3.545654703775125*^9, 3.54565470647825*^9}, {3.5461194294657474`*^9, 
   3.5461194300747824`*^9}, 3.5461238282983465`*^9, {3.546179309788824*^9, 
   3.546179310207848*^9}, 3.546254187470374*^9, 3.5462911198194304`*^9, {
   3.546336019492545*^9, 3.54633602097863*^9}, {3.5466026871642013`*^9, 
   3.546602701173026*^9}, {3.5466027390186925`*^9, 3.546602740906296*^9}, {
   3.5466030856825013`*^9, 3.5466031429190016`*^9}, 3.5466036838031516`*^9, 
   3.546603838758224*^9, {3.5466039298155837`*^9, 3.5466039303615847`*^9}, {
   3.546603989532489*^9, 3.546603992699294*^9}, {3.5466041566867824`*^9, 
   3.546604170711207*^9}, 3.546609570377092*^9, {3.5466129529466333`*^9, 
   3.546612953461434*^9}, {3.5466165615449715`*^9, 3.5466165620909724`*^9}, {
   3.546621215807147*^9, 3.546621216353148*^9}, {3.5467036651371684`*^9, 
   3.5467037049952383`*^9}, 3.546709066913888*^9, {3.54671335965131*^9, 
   3.5467133599165106`*^9}, 3.546716626294013*^9, 3.546718313560519*^9, 
   3.5467199103918524`*^9, {3.546720141385064*^9, 3.546720142650137*^9}, 
   3.546723236623102*^9, {3.5467250496548014`*^9, 3.5467250501038275`*^9}, 
   3.5467270687272863`*^9, 3.546729257230461*^9, {3.5468977714356284`*^9, 
   3.5468977840872507`*^9}}]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Integrate 2-photon Lineshape over Doppler Distribution", "Subtitle",
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"transProb", "[", 
     RowBox[{"\[Omega]1_", ",", "\[Omega]2_", ",", "v_"}], "]"}], "=", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"fDPrime", "=", "fMinFin"}], "fMaxFin"], 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"f", "=", "fMinGnd"}], "fMaxGnd"], 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          FractionBox[
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"2", " ", "\[Pi]", " ", 
                RowBox[{"\[Nu]GndToFin", "[", 
                 RowBox[{"f", ",", "fDPrime"}], "]"}]}], "-", 
               RowBox[{"\[Omega]1", 
                RowBox[{"(", 
                 RowBox[{"1", "+", 
                  FractionBox["v", "c"]}], ")"}]}], "-", 
               RowBox[{"\[Omega]2", 
                RowBox[{"(", 
                 RowBox[{"1", "-", 
                  FractionBox["v", "c"]}], ")"}]}]}], ")"}], "2"], "+", 
            FractionBox[
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"]}],
            ")"}]}], "+", "\[IndentingNewLine]", 
         RowBox[{
          FractionBox[
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"], 
          "/", 
          RowBox[{"(", 
           RowBox[{
            SuperscriptBox[
             RowBox[{"(", 
              RowBox[{
               RowBox[{"2", " ", "\[Pi]", " ", 
                RowBox[{"\[Nu]GndToFin", "[", 
                 RowBox[{"f", ",", "fDPrime"}], "]"}]}], "-", 
               RowBox[{"\[Omega]1", 
                RowBox[{"(", 
                 RowBox[{"1", "+", 
                  FractionBox["v", "c"]}], ")"}]}], "-", 
               RowBox[{"\[Omega]2", 
                RowBox[{"(", 
                 RowBox[{"1", "+", 
                  FractionBox["v", "c"]}], ")"}]}]}], ")"}], "2"], "+", 
            FractionBox[
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"\[Gamma]Fin", "+", "\[Gamma]T"}], ")"}], "2"], "4"]}],
            ")"}]}]}], ")"}], 
       RowBox[{
        UnderoverscriptBox["\[Sum]", 
         RowBox[{"m", "=", 
          RowBox[{"-", "f"}]}], "f"], 
        RowBox[{
         UnderoverscriptBox["\[Sum]", 
          RowBox[{"mDPrime", "=", 
           RowBox[{"-", "fDPrime"}]}], "fDPrime"], 
         SuperscriptBox[
          RowBox[{"Abs", "[", 
           RowBox[{
            UnderoverscriptBox["\[Sum]", 
             RowBox[{"fPrime", "=", "fMinInt"}], "fMaxInt"], 
            RowBox[{
             UnderoverscriptBox["\[Sum]", 
              RowBox[{"mPrime", "=", 
               RowBox[{"-", "fPrime"}]}], "fPrime"], 
             RowBox[{
              RowBox[{"R", "[", 
               RowBox[{"length", ",", "v", ",", "f", ",", "fPrime"}], "]"}], 
              "\[IndentingNewLine]", 
              RowBox[{"(", 
               RowBox[{
                FractionBox[
                 RowBox[{"\[Gamma]Int", "+", "\[Gamma]T"}], "2"], "/", 
                RowBox[{"(", 
                 RowBox[{
                  FractionBox[
                   RowBox[{"\[ImaginaryI]", " ", 
                    RowBox[{"(", 
                    RowBox[{"\[Gamma]Int", "+", "\[Gamma]T"}], ")"}]}], "2"], 
                  "+", 
                  RowBox[{"(", 
                   RowBox[{
                    RowBox[{"\[Omega]1", 
                    RowBox[{"(", 
                    RowBox[{"1", "+", 
                    FractionBox["v", "c"]}], ")"}]}], "-", " ", 
                    RowBox[{"2", " ", "\[Pi]", " ", 
                    RowBox[{"\[Nu]GndToInt", "[", 
                    RowBox[{"f", ",", "fPrime"}], "]"}]}]}], ")"}]}], ")"}]}],
                ")"}], 
              SubscriptBox["dGndToInt", 
               RowBox[{"f", ",", "fPrime", ",", "m", ",", "mPrime"}]], 
              SubscriptBox["dIntToFin", 
               RowBox[{
               "fPrime", ",", "fDPrime", ",", "mPrime", ",", 
                "mDPrime"}]]}]}]}], "]"}], "2"]}]}]}]}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"fnToIntegrate", "[", 
     RowBox[{"fR_", ",", "v_"}], "]"}], "=", 
    RowBox[{
     UnderoverscriptBox["\[Sum]", 
      RowBox[{"n1", "=", "1"}], 
      RowBox[{"Length", "[", "n1List", "]"}]], 
     RowBox[{
      UnderoverscriptBox["\[Sum]", 
       RowBox[{"n2", "=", "1"}], 
       RowBox[{"Length", "[", "n2List", "]"}]], 
      RowBox[{
       RowBox[{"dopplerDist", "[", "v", "]"}], 
       RowBox[{"transProb", "[", 
        RowBox[{
         RowBox[{"\[Omega]L", "[", 
          RowBox[{
           RowBox[{
           "n1List", "\[LeftDoubleBracket]", "n1", "\[RightDoubleBracket]"}], 
           ",", "fR"}], "]"}], ",", 
         RowBox[{"\[Omega]L", "[", 
          RowBox[{
           RowBox[{
           "n2List", "\[LeftDoubleBracket]", "n2", "\[RightDoubleBracket]"}], 
           ",", "fR"}], "]"}], ",", "v"}], "]"}]}]}]}]}], ";"}], 
  "\[IndentingNewLine]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vIntSpec", "[", 
    RowBox[{"fRep_", "?", "NumericQ"}], "]"}], ":=", 
   RowBox[{"NIntegrate", "[", 
    RowBox[{
     RowBox[{"fnToIntegrate", "[", 
      RowBox[{"fRep", ",", "v"}], "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"v", ",", 
       RowBox[{"-", "60000"}], ",", "60000"}], "}"}], ",", 
     RowBox[{"MinRecursion", "\[Rule]", "3"}], ",", 
     RowBox[{"MaxRecursion", "\[Rule]", "10"}], ",", 
     RowBox[{"AccuracyGoal", "\[Rule]", "4"}], ",", 
     RowBox[{"PrecisionGoal", "\[Rule]", "2"}]}], "]"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.538523548469115*^9, 3.5385235524603434`*^9}, {
   3.5385832601238136`*^9, 3.538583261230877*^9}, {3.540322004734375*^9, 
   3.54032201065625*^9}, {3.5418540524747343`*^9, 3.541854054143945*^9}, {
   3.5418542851034255`*^9, 3.5418542880986447`*^9}, {3.5456539941345*^9, 
   3.5456540111345*^9}, {3.54565470954075*^9, 3.545654710025125*^9}, {
   3.5468979843292027`*^9, 3.54689799428202*^9}, 3.5468985603758144`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fnToIntegrate", "[", 
  RowBox[{"fRepAvg", ",", "0"}], "]"}]], "Input",
 CellChangeTimes->{{3.5468985440893855`*^9, 3.5468985511717978`*^9}}],

Cell[BoxData["8.47457245323287`*^-6"], "Output",
 CellChangeTimes->{{3.5468985515149984`*^9, 3.5468985634178195`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Timing", "[", 
  RowBox[{"vIntSpec", "[", "fRepAvg", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.5461162527600503`*^9, 3.546116274539296*^9}, {
  3.5461164405267897`*^9, 3.5461164434699583`*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"29.671999999999997`", ",", "5.28446339816896`"}], "}"}]], "Output",
 CellChangeTimes->{{3.5461162499938917`*^9, 3.5461162900551834`*^9}, {
   3.5461164416948566`*^9, 3.546116450708372*^9}, 3.546120878103605*^9, 
   3.546123979689005*^9, 3.5461616669507113`*^9, 3.5461797046944113`*^9, 
   3.54625434844304*^9, 3.546291275643343*^9, 3.5463361787606544`*^9, 
   3.5468980169800596`*^9, 3.546898594103073*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Calculation", "Subtitle",
 CellChangeTimes->{{3.456681446140625*^9, 3.45668144753125*^9}, {
  3.5379007941050625`*^9, 3.537900798020669*^9}},
 FontColor->RGBColor[0, 0.356863, 0]],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Delta]f", "=", "0.25"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fr", "=", "fRepMin"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"fRepList", "=", 
    RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Form", " ", "the", " ", "fRepList"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"While", " ", "[", 
  RowBox[{
   RowBox[{"fr", "<", "fRepMax"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"fRepList", "=", 
      RowBox[{"Append", "[", 
       RowBox[{"fRepList", ",", "fr"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"fr", "+=", "\[Delta]f"}], ";"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{
  3.456746939*^9, {3.5372003913424997`*^9, 3.5372003937643747`*^9}, {
   3.5378885919057503`*^9, 3.5378885935932503`*^9}, {3.537901719773398*^9, 
   3.53790172663741*^9}, 3.5379017983975363`*^9, 3.5385235615778646`*^9, 
   3.538524008413422*^9, {3.5401635046948967`*^9, 3.5401635052096977`*^9}, 
   3.540322016265625*^9, {3.54185313798652*^9, 3.5418531382305336`*^9}, 
   3.5461162162409616`*^9, {3.546898926275457*^9, 3.5468989270086584`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"fRepMax", "-", "fRepMin"}]], "Input",
 CellChangeTimes->{{3.546116314559585*^9, 3.5461163148105993`*^9}}],

Cell[BoxData["1500.`"], "Output",
 CellChangeTimes->{{3.546116310488352*^9, 3.5461163151786203`*^9}, 
   3.5461164507423744`*^9, 3.546120878143607*^9, 3.5461239797210073`*^9, 
   3.5461616670037146`*^9, 3.546179704762415*^9, 3.5462543484586396`*^9, 
   3.546291275697346*^9, 3.546336178839659*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"fRepList", ";"}]], "Input",
 CellChangeTimes->{3.5461613531097608`*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Export", "[", 
     RowBox[{
     "\"\<C:\\\\RbDFCS\\\\Rb87_5D52.txt\>\"", ",", "totalSpec", ",", 
      "\"\<Table\>\""}], "]"}], ";"}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"tt", "=", 
     RowBox[{"AbsoluteTime", "[", "]"}]}], ";"}], "\[IndentingNewLine]", " ", 
   RowBox[{"(*", 
    RowBox[{
    "distribute", " ", "NIntegrate", " ", "among", " ", "CPU", " ", "cores"}],
     "*)"}], "\n", 
   RowBox[{
    RowBox[{"specList", "=", 
     RowBox[{"ParallelTry", "[", 
      RowBox[{"vIntSpec", ",", "fRepList", ",", 
       RowBox[{"Length", "[", "fRepList", "]"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"ttt", "=", 
     RowBox[{
      RowBox[{"AbsoluteTime", "[", "]"}], "-", "tt"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"Print", "[", 
    RowBox[{"\"\<Execution time = \>\"", ",", "ttt"}], "]"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"fRepList", "=", "fRepList"}], ";"}], " ", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"totalSpec", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         FractionBox[
          RowBox[{
          "fRepList", "\[LeftDoubleBracket]", "t", "\[RightDoubleBracket]"}], 
          SuperscriptBox["10", "6"]], ",", 
         RowBox[{"Chop", "[", 
          RowBox[{
          "specList", "\[LeftDoubleBracket]", "t", "\[RightDoubleBracket]"}], 
          "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"t", ",", "1", ",", 
         RowBox[{"Length", "[", "specList", "]"}]}], "}"}]}], "]"}]}], ";"}], 
   "\[IndentingNewLine]", 
   RowBox[{"ListPlot", "[", 
    RowBox[{"totalSpec", ",", " ", 
     RowBox[{"PlotRange", "\[Rule]", " ", "All"}], ",", " ", 
     RowBox[{"Joined", "\[Rule]", " ", "True"}]}], " ", "]"}]}]}]], "Input",
 CellChangeTimes->{{3.5385255757270675`*^9, 3.5385255896838655`*^9}, 
   3.538525623446797*^9, {3.5385833530161266`*^9, 3.5385833644487805`*^9}, 
   3.5386070295553474`*^9, {3.5391281630166893`*^9, 3.539128167970973*^9}, {
   3.54006852448602*^9, 3.5400685308040204`*^9}, {3.5410239000592422`*^9, 
   3.5410239015417175`*^9}, {3.5410239950156674`*^9, 
   3.5410239980430374`*^9}, {3.5410240329982376`*^9, 3.541024060104122*^9}, {
   3.5410262133987827`*^9, 3.541026267538763*^9}, {3.5461220657525344`*^9, 
   3.546122073575982*^9}}],

Cell[BoxData[
 InterpretationBox[
  RowBox[{"\<\"Execution time = \"\>", "\[InvisibleSpace]", 
   "14006.0060977`11.597859304583833"}],
  SequenceForm["Execution time = ", 14006.0060977`11.597859304583833],
  Editable->False]], "Print",
 CellChangeTimes->{3.5420852997627926`*^9, 3.5424503050821247`*^9, 
  3.545654925150125*^9, 3.545685493697*^9, 3.546116180320907*^9, 
  3.546117521028591*^9, 3.546121956646294*^9, 3.5461250719814806`*^9, 
  3.5461763286113105`*^9, 3.546192811402072*^9, 3.54626686352605*^9, 
  3.546305089013423*^9, 3.546350184863758*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"Export", "[", 
   RowBox[{"oFile", ",", "totalSpec", ",", "\"\<Table\>\""}], "]"}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.5385255757270675`*^9, 3.5385255896838655`*^9}, 
   3.538525623446797*^9, {3.5385833530161266`*^9, 3.5385833644487805`*^9}, 
   3.5386070295553474`*^9, {3.5391281630166893`*^9, 3.539128167970973*^9}, {
   3.5401635420881624`*^9, 3.5401635432581644`*^9}, {3.54023048903125*^9, 
   3.54023049575*^9}, {3.5418531660981274`*^9, 3.5418531666921616`*^9}}]
}, Open  ]]
},
Evaluator->"Local",
WindowSize->{1132, 717},
WindowMargins->{{64, Automatic}, {27, Automatic}},
Magnification->1.350000023841858,
FrontEndVersion->"7.0 for Microsoft Windows (64-bit) (February 18, 2009)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[545, 20, 488, 13, 85, "Input"],
Cell[CellGroupData[{
Cell[1058, 37, 76, 1, 65, "Subtitle"],
Cell[1137, 40, 890, 17, 108, "Input"],
Cell[CellGroupData[{
Cell[2052, 61, 37, 0, 38, "Subsubtitle"],
Cell[2092, 63, 148, 5, 38, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2277, 73, 40, 0, 38, "Subsubtitle"],
Cell[2320, 75, 223, 8, 41, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[2580, 88, 107, 1, 38, "Subsubtitle"],
Cell[2690, 91, 330, 10, 41, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[3069, 107, 75, 1, 65, "Subtitle"],
Cell[3147, 110, 102, 1, 65, "Subtitle"],
Cell[3252, 113, 149, 2, 38, "Subsubtitle"],
Cell[3404, 117, 271, 6, 57, "Input"],
Cell[CellGroupData[{
Cell[3700, 127, 105, 1, 38, "Subsubtitle"],
Cell[3808, 130, 594, 12, 58, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4439, 147, 101, 1, 38, "Subsubtitle"],
Cell[4543, 150, 219, 6, 38, "Input"],
Cell[4765, 158, 167, 5, 38, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[4981, 169, 106, 1, 65, "Subtitle"],
Cell[5090, 172, 222, 3, 38, "Subsubtitle"],
Cell[5315, 177, 225, 5, 57, "Input"],
Cell[CellGroupData[{
Cell[5565, 186, 199, 3, 38, "Subsubtitle"],
Cell[5767, 191, 204, 5, 57, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6008, 201, 109, 1, 38, "Subsubtitle"],
Cell[6120, 204, 301, 8, 62, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6458, 217, 389, 5, 38, "Subsubtitle"],
Cell[6850, 224, 948, 17, 38, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7835, 246, 413, 6, 38, "Subsubtitle"],
Cell[8251, 254, 554, 9, 38, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[8854, 269, 115, 1, 65, "Subtitle"],
Cell[8972, 272, 249, 3, 38, "Subsubtitle"],
Cell[9224, 277, 250, 6, 57, "Input"],
Cell[CellGroupData[{
Cell[9499, 287, 199, 3, 38, "Subsubtitle"],
Cell[9701, 292, 275, 6, 57, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10013, 303, 115, 1, 38, "Subsubtitle"],
Cell[10131, 306, 347, 9, 62, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10515, 320, 252, 3, 38, "Subsubtitle"],
Cell[10770, 325, 630, 16, 58, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11437, 346, 315, 4, 38, "Subsubtitle"],
Cell[11755, 352, 1232, 22, 62, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[13024, 379, 350, 5, 38, "Subsubtitle"],
Cell[13377, 386, 1048, 19, 38, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14462, 410, 371, 5, 38, "Subsubtitle"],
Cell[14836, 417, 1040, 19, 38, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[15925, 442, 155, 2, 65, "Subtitle"],
Cell[16083, 446, 240, 3, 38, "Subsubtitle"],
Cell[16326, 451, 323, 7, 57, "Input"],
Cell[CellGroupData[{
Cell[16674, 462, 199, 3, 38, "Subsubtitle"],
Cell[16876, 467, 370, 7, 57, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17283, 479, 104, 1, 38, "Subsubtitle"],
Cell[17390, 482, 347, 9, 62, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17774, 496, 293, 4, 38, "Subsubtitle"],
Cell[18070, 502, 764, 16, 58, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[18883, 524, 139, 1, 65, "Subtitle"],
Cell[19025, 527, 175, 2, 38, "Subsubtitle"],
Cell[19203, 531, 442, 15, 38, "Input"],
Cell[19648, 548, 1602, 47, 114, "Input"],
Cell[21253, 597, 2287, 55, 292, "Input"],
Cell[CellGroupData[{
Cell[23565, 656, 57, 0, 38, "Subsubtitle"],
Cell[23625, 658, 1384, 45, 98, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25046, 708, 97, 1, 38, "Subsubtitle"],
Cell[25146, 711, 1254, 33, 174, "Input"],
Cell[26403, 746, 1184, 32, 147, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27624, 783, 115, 1, 38, "Subsubtitle"],
Cell[27742, 786, 2134, 54, 374, "Input"],
Cell[29879, 842, 2166, 55, 397, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[32094, 903, 89, 1, 65, "Subtitle"],
Cell[32186, 906, 521, 12, 65, "Input"],
Cell[CellGroupData[{
Cell[32732, 922, 280, 7, 57, "Input"],
Cell[33015, 931, 423, 6, 38, "Output"]
}, Open  ]],
Cell[33453, 940, 222, 7, 62, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33712, 952, 73, 1, 65, "Subtitle"],
Cell[33788, 955, 39, 0, 38, "Subsubtitle"],
Cell[33830, 957, 681, 13, 41, "Input"],
Cell[CellGroupData[{
Cell[34536, 974, 102, 1, 38, "Subsubtitle"],
Cell[34641, 977, 2116, 36, 67, "Input"],
Cell[36760, 1015, 198, 5, 57, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[36995, 1025, 90, 1, 38, "Subsubtitle"],
Cell[37088, 1028, 507, 15, 107, "Input"],
Cell[37598, 1045, 833, 25, 107, "Input"],
Cell[CellGroupData[{
Cell[38456, 1074, 1166, 26, 38, "Input"],
Cell[39625, 1102, 813, 14, 62, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[40475, 1121, 1215, 26, 38, "Input"],
Cell[41693, 1149, 644, 11, 38, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[42386, 1166, 102, 1, 38, "Subsubtitle"],
Cell[42491, 1169, 357, 10, 38, "Input"],
Cell[CellGroupData[{
Cell[42873, 1183, 101, 1, 95, "Section"],
Cell[42977, 1186, 382, 9, 41, "Input"],
Cell[43362, 1197, 416, 12, 41, "Input"],
Cell[43781, 1211, 736, 16, 58, "Input"],
Cell[CellGroupData[{
Cell[44542, 1231, 1302, 36, 103, "Input"],
Cell[45847, 1269, 874, 18, 43, "Output"]
}, Open  ]],
Cell[46736, 1290, 514, 12, 64, "Input"],
Cell[47253, 1304, 785, 22, 89, "Input"],
Cell[CellGroupData[{
Cell[48063, 1330, 379, 12, 38, "Input"],
Cell[48445, 1344, 957, 15, 38, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[49439, 1364, 284, 7, 62, "Input"],
Cell[49726, 1373, 3036, 53, 313, "Output"]
}, Open  ]],
Cell[52777, 1429, 2112, 44, 85, "Input"]
}, Open  ]]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[54950, 1480, 112, 1, 65, "Subtitle"],
Cell[55065, 1483, 6383, 166, 553, "Input"],
Cell[CellGroupData[{
Cell[61473, 1653, 165, 3, 38, "Input"],
Cell[61641, 1658, 118, 1, 40, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[61796, 1664, 221, 4, 38, "Input"],
Cell[62020, 1670, 446, 7, 38, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[62515, 1683, 185, 3, 65, "Subtitle"],
Cell[62703, 1688, 1192, 27, 177, "Input"],
Cell[CellGroupData[{
Cell[63920, 1719, 129, 2, 38, "Input"],
Cell[64052, 1723, 297, 4, 38, "Output"]
}, Open  ]],
Cell[64364, 1730, 95, 2, 38, "Input"],
Cell[CellGroupData[{
Cell[64484, 1736, 2411, 61, 244, "Input"],
Cell[66898, 1799, 558, 10, 28, "Print"]
}, Open  ]],
Cell[67471, 1812, 509, 9, 38, "Input"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
